Mesoscopic quantum switching of a Bose-Einstein condensate in an optical lattice 
governed by the parity of the number of atoms 



On 
O 
O 
(N 



oo 
(N 



V. S. Shchesnovich 

Centra de Ciencias Naturais e Humanas, Universidade Federal do ABC, Santo Andre, SP, 09210-170 Brazil 

It is shown that for a iV-boson system the parity of N can be responsible for a qualitative difference 
in the system response to variation of a parameter. The nonlinear boson model is considered, which 
describes tunneling of boson pairs between two distinct modes Xi : 2 of the same energy and applies to 
a Bose-Einstein condensate in an optical lattice. By varying the lattice depth one induces the parity- 
dependent quantum switching, i.e. X\ — > X2 for even N and X\ — » Xi for odd N, for arbitrarily 
large N. A simple scheme is proposed for observation of the parity effect on the mesoscopic scale 
by using the bounce switching regime, which is insensitive to the initial state preparation (as long 
as only one of the two Xi modes is significantly populated), stable under small perturbations and 
requires an experimentally accessible coherence time. 
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Mesoscopic quantum phenomena lie in between the big 
and small: the macroscopic classical world and the mi- 
croscopic quantum world. The Bose-Einstein condensate 
(BEC) is such a mesoscopic effect, i.e. a big "matter- 
wave" . An order parameter governed by the Gross- 
Pitaevskii equation [l|, is usually attributed to BEC. 
The order parameter corresponds to the mean-field the- 
ory, i.e. to the limit of large number of bosons: N — > 00 
at a constant density. The latter, on the other hand, 
is equivalent to the classical limit of the discrete WKB 
approach [J], with l/N playing the role of the Planck 
constant (see, for instance, Refs. @,@|). 

The mean-field limit is a singular limit of the full 
quantum description and suffers from deficiency, e.g. at 
a dynamic instability due to the back-reaction of 
the quantum fluctuations ML or the appearance of the 
Schrodinger cat-like states |8j. In this connection one 
can mention the "even-odd" effect, first predicted for the 
spin systems [§] and observed in the small (S ~ 10) mag- 
netic molecular clusters as the parity-dependent tunnel- 
ing splitting [HJ. The parity effect was also found in the 
decay of the Josephson 7T-states and in the boson- 
Josephson model [12j. The tunneling splitting, however, 
decreases exponentially in N, for N ^> 1, restricting its 
observation to the sub-mesoscopic scale. One may won- 
der whether it is possible to magnify the microscopic par- 
ity difference to a mesoscopic scale and how? Such an 
effect would be also an interesting manifestation of the 
singularity in the N — > 00 limit of the discrete WKB. 

The aim of this rapid communication is to present a 
solution: one must look for a dynamic parity effect which 
allows for a massive constructive quantum interference. 
Moreover, the feasibility of the experimental observation 
is shown. The mesoscopic parity effect appears in the 
response to variation of a parameter in the nonlinear 
two- mode boson model [H, a nonlinear variant of the 
celebrated boson- Josephson model HI, 3j- 

The nonlinear two-mode boson model is formulated as 
follows. Suppose that a single-particle Hamiltonian Hq 



has two equal energy states X\ and Xi and that the inter- 
action term H mi = § J d 3 x^t (x^ (x)ip(x)t/j(x) in the 
many-body Hamiltonian H = j d 3 x-0t ( x )iJ o -0(x) + H lnt 
is smaller than the energy gap of H isolating the reso- 
nant subspace. Projecting on the resonant states, tp{x.) = 
b iVXi ( x ) + b 2 i Px 2 ( x )i °ne arrives at the Hamiltonian: 

H m t = § {j2Xti,t2,n,j2 b l 1 b l 2 b j 1 b j 2 }> e {1,2}, 

where Xii,i 2 ,jx,j2 = Jd^^^^^- We con- 
sider the situation when the bosons hop between the 
modes Xi by pairs, i.e. when Xi2jj — f° r 3 = 
This type of coupling describes the intraband tunneling 
of BEC in a square optical lattice where the reso- 
nant modes are the high-symmetry points of the Brillouin 
zone, X\ — (fcs,0) and X2 = (0, fcs), and the quasimo- 
mentum conservation makes Xi2jj vanish. Moreover, it 
also applies to BEC on a rotated ring lattice [l5j- The 
intraband BEC tunneling Hamiltonian reads 
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where rij = and A = xma/Xim (0 < A < 1) is 
the only parameter in the model (see for details Ref. 
14j|). The corresponding Schrodinger equation is cast as 
ihd T \%) = H\V), where h = 2/N is the effective Planck 
constant and the dimensionless time r = {2gNxim/K)t, 
with g = Airh 2 a s /m, depends on N through the density 
only. 



Hamiltonian (fTJ) features 14| a quantum phase transi- 
tion at the top of the spectrum related to the mean-field 
symmetry-breaking bifurcation between the stationary 
point ((b\bi)/N = 1/2,0 = arg((4) 2 &?) = 0), corre- 
sponding to the equally populated X\ modes, which is 
stable for A > A c = 1/3, and the selftrapping stationary 
points {{b\bi)/N < 1 or {b\b 2 )/N < 1, <f> undefined), 
stable for A < A c . It also has a parity-dependent en- 
ergy spectrum (see also Fig Q3a)). There are two in- 
variant subspaces corresponding to the even and odd 
occupation numbers k in the Fock basis, i.e. |^/) = 
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y,*- n Ck\k,N-k), with \k,N-k) = ^MHJllvac). 

The projections of H on the even (s = or "cv") and 
odd (s = 1 or "od") subspace are given as 
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where, respectively for even and odd N: L = 
N/2 and L = (N - l)/2 in the case of H^ cv \ 
while L = N/2 and L = (N + l)/2 in the case 
of i/(° d '. Here a fe = (2A - 1) A (!_.&) and 

fa = f [# (# + 7f) i 1 ~ #) C 1 _ # + £)] 1/2 ' 

When AT is odd we have PH^P = #( od \ where 
P = diag(l, . . . , 1) T (i.e. the transposed 1). In this case 
P\Ef v) ) = |^ (od) ), i.e. the energy levels are doubly de- 
generate. On the other hand, this is a consequence of the 
Kramers theorem [l6j . Indeed, Hamiltonian |T]) is equiv- 
alent to a spin model H s = (1 - A)Sf + 2AS^ with the 
total spin S = N/2, if we associate S x = (b\b2 + b\b\)/2, 
S y = {b\b 2 - b%)/2i and S z = (b{bi - blh)/2. 

When A" is even the projected Hamiltonian is in- 
variant under the exchange symmetry, i.e. PH^P = 
H( s \ s € {0, 1}. One control parameter can not cause 
the energy level degeneracy [17] . thus the eigenvectors 
of must satisfy the exchange symmetry, namely 
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A > the eigenvalues of appear in the 
form of very narrow doublets (see also Fig. QIa)) due to 
the "selftrapping states" being strongly localized at the 
respective Xi mode (e.g. typically | (fc, A^ — fc|^) | 2 < 10~ 6 
for k > 25 14]). Consequently, each also has the 
quasi-degenerate spectra (since ij( ov )— ij( od ) oc A) which 
become finer for larger deviations A c — A, exactly as the 
numerics indicates. 

Consider now the following experimentally realizable 
setup: initially just one of the Xi modes is significantly 
populated (to achieve this one can use the non-adiabatic 
loading [l8[ into one of the two resonant Bloch states of 
the lattice with Ai < 1/3). By varying the lattice pa- 
rameter (e.g. by changing the lattice depth) between Ai 
and A2, with A2 > A c , one drives the system across the 
phase transition and back to force a switching-like dy- 
namics between the selftrapping states at the Xi modes. 
Remarkably, for a general initial state localized at just 
one Xi mode, e.g. = J2k^> N Ck\k,N—k), where the 
distribution Ck is not important, the result qualitatively 
depends on the parity of AT, see Fig. [ljb). Note that the 
switching is between the Bloch modes with orthogonal 
Bloch vectors: X\ = (fee, 0) and X 2 — (0, fee). 

In the adiabatic limit the mechanism of the switch- 
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FIG. 1: (Color online) Panel (a): the energy level structure 
of H {cv) for N = 40 (solid lines) and 41 (dashed lines). Panel 
(b): the average ratio of the Xi-mode occupation number 
(as) = (b}bi)/N (dark solid lines) and (as) - ^, (x) + 
(light dashed lines), with Ax = {{btbx/N - (x}) 2 } 1/2 . The 
upper lines (solid and dashed) correspond to N — 40 and 
the lower ones to N = 41. The initial state is a Gaus- 
sian |*) = A (T J]f =0 e- fe2/CT2 |A;,7V-fc>, where a = 4. Here 
A(t) = Ai + 1/2(A 2 - Ai) [tanh(s[r - n]) - tanh («[r - r 2 ])] 
with k = 10~ 4 , n = 7000 and r 2 = 12800. 




FIG. 2: (Color online) The probabilities {E^ |0, N) | 2 , 
j = L + 1 and j — L, given for N — 40 by the solid 
and dashed lines (almost coincide), for N = 60 by the 
"+" and "x", for N = 80 by the "□" and "o" and 
for N — 100 by the "*" and "o" symbols. The upper 
dash-dotted line gives Y? m =a \{E££{_ m \0, N)\ 2 . The upper 
and lower dotted lines give, respectively, \(E^\\0, N}\ 2 and 

E™=i K-^iTi-ml !-^)! 2 for odd N - The data > connected by 
the dotted lines to guide the eye, that lie off the central curve 
are due to the quasi-degeneracy on the order of round-off er- 
ror. 



ing for a localized initial distribution Ck can be under- 
stood by considering separately the even and odd invari- 
ant subspaces. For simplicity, consider the initial state 
|*) = |0, N) £ H^l (the case of |W) = |1, N-l) € ff (od) 
is similar). Fig. [2] shows that there are but few sig- 
nificant terms (from the top of the spectrum) in the 
expansion of the initial state over the eigenvectors (re- 
call the strong localization of the eigenvectors at k = 
and k — N for A < A c ). For even N there are the 
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quantum beats, i.e. oscillations of the populations, be- 
tween the successive pairs of eigenstates from the top 
of the energy spectrum. This can be understood as 
tunneling between the Xi modes. The switching oc- 
curs for the odd-7r phase differences between the first 
few pairs of the successive top energy eigenstates, i.e. 

A0 i+1 _ 2j = i Jdt [4 S li- 2i (*)-4- 2j -(*)] » (Zmj-IK 
o 

j = 0, 1, ... . For odd TV, on the other hand, the eigen- 
vectors break the P-symmetry and only those localized 
at X\ acquire nonzero amplitudes (i.e. no tunneling). 
Thus, a high visibility parity effect requires at least few 
top dynamic phase differences in both the even and odd 
subspaces (for even TV) to be odd in n, which, in the 
general case, can not be satisfied by adjusting t 2 — t\ in 

Fig. in 

The Landau-Zener-Majorana (LZM) transi- 
tions between the instantaneous energy levels 
occur for the non-adiabatic variation of A. In 
the instantaneous eigenvector basis, |^w(r)) = 



EjCj(r) exp 



o 
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Ey[r)), we have 
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dr dr 
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where it was used that (E^^E^) = {H^ does not 
have energy degeneracy) . Due to the exchange symmetry 
P, in the even TV case the LZM tunneling occurs only 

(s) 

between the levels P) with the same parity of j, whereas 
for odd TV there is a small coupling also between the 
r« _ H-00 oc i The LZM result 



adjacent levels, since i^i+i - ^2L 



N 



[lit states that |q| 2 ~ exp 



, for I 7^ j. Hence, 
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the lower limit on the adiabatic time scale, i.e. |q| 2 <C 
1, can be determined from the difference between the 



top energy levels E^ - En 1 cx jfa for A -> A c [l4j 
(see also below). We get r a d 3> TV (the time scale for a 
finite phase difference A0j is r p h ~ TV, since h = 2/N). 
Therefore, the adiabatic case also requires an extremely 
long coherence time and thus it is not realistic at all. 

In the other limit, Fig. [31 when A(r) is a step-like 
function (e.g. similar to the one used in Fig. [ljb), but 
with k ~ 1 or larger) taking two values Ai < A c and 
A2 = A c , the switching regime, the bounce switching for 
below, has all the needed properties. 

Indeed, the bounce switching possesses the TV-scaling 
property, thus it survives in the N — > 00 limit, and is 
insensitive to the initial distribution Ck (localized at just 
one X-mode) . These features originate from the fact that 
all the dynamic phases are odd in tt. Indeed, Hamil- 
tonian {TJ) in the coherent basis ai,2 = (61 T ib 2 )/\/2 
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FIG. 3: (Color online) The bounce switching for the initial 

state |*) = A„Y<k=o e ~ k2/ ' y2 \ k > N ~ k )- Panel ( a ) N = 200 
and panel (b) TV = 201. In both panels, the average ratio (x) 
(the thick lines) with the side lines (a;) — ^ and (x) + ^ 
(the thin lines) are given for a = 0.057V (dash-dotted lines), 
a — 0.17V (dashed lines), and a = 0.27V" (solid lines). Here 
A(r) is as in Fig. [TJ with k = 8, Ai = 0.25, A 2 = 1/3, and 

T'2 — Tl = 37V7T. 



I = 1, [y] + 1. Hence, the dynamic phase differences 

are given as A6^ +1 _ s _ 2n = ^Ar c , n = 0, 1, 2, . . ., 
where Ar c is the hold time at A c . By setting At c = 37rTV 
one obtains A9^ +1 _ s _ 2n = {An + 2s + l)ir. 
The actual physical time, t = t p h-^ with t p h ~ 



. c j -igxxxxi ' 

is TV-independent due to the scaling property Ar c = 
37r7V. For a condensate in a square lattice of n cites 

of size d in the tight transverse trap with the oscilla- 

.2 

tor length a± the coherence time is At ~ h ^ ±n ■ For 
87 Rb, n = 64, d — 0.5/im and a± = 0.1/xm the re- 
quired coherence time is At ~ 0.3 s. The switching can 
be detected by releasing the optical lattice and observ- 
ing the direction of the interference pattern. For the lat- 
tice V = Vq [cos(2k B x) + cos(2k B y)] the parameter range 
is V = 0.58P fl for A = 0.25 and V = 0M73E R for 
A = 1/3. Finally, the applicability condition for the two- 
mode model Pnl/P/?. -C 1, where E R = h 2 k B /2m and 
Pnl = ffA^xim, can be cast as TV -C na±/a s . 

Experimental observation requires stability of the dy- 
namic parity effect under perturbations. To have an idea 
on the bounce switching stability, consider first the gen- 
eral perturbation within the two-mode model: 

ffpert = ^(n a - n 2 ) + — {b\b 2 + b\b x ), (4) 
where e and J (given in terms of the nonlinear energy 
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-Enl) account for the imperfections of the optical lattice 
and for a magnetic trap (but, see below). The nonlinear 
interaction terms gxi2jj{b\b\[b\ + &|] + h.c.) discarded 
when deriving Eq. JT]) also reduce to the J-part in Eq. 
(g|) with J = xiw/XUU- The first term in Eq. @ 
preserves the decomposition of the model as in Eq. |2|), 
whereas the second one breaks it and both break the 
exchange symmetry P. In the a-operator basis H pcrt = 
jj(a\a2+a\ai) — ^(a\a,2 — a\a{) and the matrix elements 
of -f/pert between the eigenstates of H are much smaller 
than the energy differences at A c if \e — iJ\ <C 1/N 2 . 
However, extensive numerical simulations show that for 
N < (e 2 + J 2 )~ 1//4 the system still exhibits essentially the 
same parity effect as in Fig. H For N < |A 2 - 1/3|~ 1/2 , 
the parity effect is also insensitive to an imprecise tuning 
of A2 in Fig. [3] to the critical value. 

An additional weak magnetic trap Vt Y = mu 2 (i 2 + 
y 2 )/2 + V± (z) is always a part of the experimental setup, 
leading to the J-term in Eq. (|4|). However, this con- 
tribution is exponentially small. Indeed, for a weak 



trap 



h/muj ~ nd 3> d the unperturbed Bloch 



wave <Pk(x) acquires a factor given by a product of a 
polynomial and a Gaussian in (x,y), i.e. <f\^ (x) = 
G(z)P(x, y) exp{— x ^ }ipk(x). The Gaussian factor 
defines the order of the non-diagonal matrix element of 
Vt r between the two resonant Bloch waves. We have 



d Xt ^(fe B ,0) l/ tr^ (0i/! . £ 



hm exp 



Ad 2 



(5) 



therefore J in Eq. dU reads J ~ (Huj/Enl) exp{— 7r 2 n/4}. 
To bound the pre-exponential factor hw/E^ ~ 
nd 2 a±/(£ 2 a s N) ~ a±/(a s N) one needs JV ~ a±/a s com- 
patible with the applicability condition N <C na±/a s . 

There still remains to consider the transitions to the 
non-resonant modes, due to the nonlinear term of the 
full boson Hamiltonian. The latter, however, preserve 
the quasi-momcntum (with the exponentially small cor- 
rection due to the magnetic trap) and, hence, the par- 
ity of N, since the bosons leave the resonant modes by 
pairs. More detrimental than the setup imperfections 
considered above is the loss of atoms, for instance the 
scattering of BEC atoms with the cloud of hot atoms, 
which will wash out the parity effect. To prevent this, a 
smaller then in a usual BEC number of cold atoms can 
be used, e.g. N on the order of few hundred atoms. 

In conclusion, for a flexible control parameter, the non- 
linear two-mode boson model possesses a bounce switch- 
ing regime with the qualitatively different outcome of 
switching for even and odd N. This regime is insensitive 
to the initial state preparation (with just one resonant 
mode being significantly populated), shows stability to 
small perturbations and requires an experimentally ac- 
cessible coherence time, thus allowing for observation of 



the even-odd effect on the mesoscopic scale. As a general 
perspective, one can observe that the nonlinear two-mode 
boson model is a nonlinear variant of the two-site Bose- 
Hubbard Hamiltonian and that the second-order tunnel- 
ing applies also to the dynamics of the repulsively bound 
atom pairs in an optical lattice 2^, 21 1 and to the case of 
strong interactions reaching the fermionization limit [22| . 
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